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ABSTRACT 

In the present paper, some observations are made on the application of a 
three-dimensional viscous marching analysis (PEPSIS) to solve the glancing 
oblique shock wave/turbulent boundary layer interacticn In three-dimensions. 
The experimental data of Oskam. Vas and Bogdonoff provides a description of 
the measured flow field and is used as the basis for an extensive benchmark 
verification of the 3d visicous marching analysis. 

This paper delves Into the solution methodology, or calculatlonal proto- 
col, to examine such areas as the sensitivity of Initial conditions, mesh 
resolution, wall functions, sublayer resolution and turbulence model on the 
calculated results. 

INTRODUCTION 

The design of supersonic Inlets Is a difficult task In view of the wide 
operating range over which good performance Is desired. Design of such an 
Inlet system Is strongly affected by the requirement that the aircraft operate 
at all speeds from zero to the supersonic design point. In addition, most 
supersonic aircraft must have a supersonic cruise capability where fuel 
economy Is Important. As a result of these varied and sometimes conflicting 
requ1ren«nts, the design of a supersonic Inlet Is a difficult comprise. 

The foregoing considerations suggest that Inlet design technology would 
benefit from a detail and accurate flow field prediction procedure that In- 
cludes shock-boundary layer Interaction effects. Within rectangular Inlets 
operating at supersonic speeds, these viscous Interaction effects can be 
classified Into three types: (1) Incident shock reflections on the ramp or 

cowl, (2) glancing sidewall shock wave Interactions, (3) corner shock wave 
Interactions. In view of the Importance of these Interactions there have been 
a number of efforts devoted to developing. Improving and Implementing numeri- 
cal techniques to obtain solutions to the phenomena. These solutions tech- 
niques generally fall Into three categories (1) control volume methods, 

(refs. 1 to 3), (2) viscous marching methods, (refs. 4 to 6), and (3) 
Navler-Stokes solutions, (refs. 7 and 8). 

In the present paper, some observations are made on the application of the 
viscous marching analysis (ref. 6) to the glancing oblique shock wave turbu- 
lent boundary layer Interaction. An experimental study of the glancing side- 
wall Interaction was done by Oskam, Vas and Bogdonoff (ref. 9) and provides a 
detailed flow description from which to verify this analysis. It Is this 
experimental test case that forms the basis for the present study. 
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RESULTS AND DISCUSSION 

In turbulent flows, the finite difference mesh must be packed very densely 
near the wall to resolve the steep velocity gradients exsisting there. This 
heavy packing can be partially avoided by making use of the experimentally 
observed properties of turbulent boundary layers. A typical unseparated boun- 
dary layer can be divided into an inner region where the total shear stress is 
essentially constant, and an outer regTon where the total shear stress de- 
creases with distance from the wall . The inner region can be further divided 
into a laminar sublayer where the laminar shear stress dominates, a transition 
region or buffer region where both laminar and turbulent shear are in^ortant ^ 

and a fully turbulent region. * 

When wall function boundary conditions are used, the laminar sublayer is 
not resolved. Rather it is assumed that the velocity near the wall obeys the 
1 aw-of-the-wal 1 and a slip or wall function boundary condition can be de- 
fined. Using wall function boundary conditions cun in principle save substan- 
tial numbers of mesh points, particularly in three-dimensional flows since the 
near wall region is not resolved. However, in highly complex three- 
dimensional flows such as the glancing sidewall shock boundary layer inter- 
action, it is not clear if this approximation will result in proper flow 
physics being displayed. This study, therefore, delves into the use of wail 
functions as well as the companion problem of the degree of mesh lesolution 
necessary to resolve the flow physics of Interest when no-slip boundary condi- 
tions are applied. Three test cases were run with PEPSIS to examine these 
issues. A basel ine or Wall Function test case, a Coarse Mesh case with the 
same three-dimensional mesh as the Wall Function case but with no-slip 
boundary conditions and a Medium Mesh test case where the region near the 
tunnel sidewall was resolveS^ K fourth or Tunnel case was also run to demon- 
strate the complexity of the flow when the sidewall shock induced flow field 
impinges on the tunnel celling. 

PEPSIS Computer Analysis 

The PEPSIS computer code, an acronym for Parabolic Elliptic Streamwise 
Implicit Supersonic, was developed jointly by NASA Lewis Research Center and 
Scientific Research Associates. It is a fully three-dimensional computer 
design code for supersonic inlets. The computer analysis possess the general 
three-dimensional viscous nature of the Navier-Stokes equations, but takes 
advantage of realistic physical approximations to limit computer run time and 
storage requirements associated with the solution of the complete three- 
dimensional Navier-Stokes equations. The assumption is made that a "primary" i 

flow direction exists and that diffusion arising from the rate of change of 
the flow in this "primary direction" can be neglected. In this manner a set 
of steady state equations is produced for entirely supersonic flows which are • 

solved by an efficient spatial marching procedure. 

Glancing Sidewall Shock Wave Interaction 

A schematic diagram depicting the analytical and experimental test con- 
figuration used to study the glancing sidewall shock wave boundary layer 
interaction is shown in figure 1. This interaction arises when the oblique 
shock wave formed by the wedge interacts with the boundary layer exsisting 
along the tunnel wall. This interaction is fully three-dimensional since the 
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position of the shock on the tunnel M&11 surface chan 9 ts In the cross stream 
direction. Because of the skewing of the Incident shock on the tunnel wall, a 
strong transverse gradient produces a cross flow with the tunnel wall boundary 
layer in a direction along the shock wave. Downstream of the shock wave, the 
transverse static pressure gradient disappears. At a given distance normal to 
the tunnel sidewall within the boundary layer Interaction zone, the flow 
properties vary rapidly In the direction normal to the tunnel sidewall. While 
the region of study In the experimental Investigation was confined primarily 
to the center section of the tunnel sidewall, the analytical study Includes 
the ramp surface, the corner region and the tunnel sidewall. 

Figure 2 shows the shock generator along with the experimental survey 
planes, which were the same in the analysis for the purpose of data com- 
parisons. The stremnwise coordinate X was measured from the leading edge of 
the wedge, the normal coordinate Y was measured relative to the wedge surface 
and the coordinate Z was normal to the X-Y plane. For this study, the wedge 
angle was set at 10.0*. 


Calculational Test Conditions 

The coordinate system and represented computational mesh used in this 
study are depicted in figure 3 along with the identification of important 
surfaces critical to this calculation. These important surfaces Include the 
initial data plane, located upsteam of the ramp tip, the ramp surface and the 
tunnel sidewall. Also shown on figure 3 is the inviscid location of the shock 
generated by the 10.0* wedge. The analysis of the glancing sidewall shock 
wave boundary layer interaction was started upstream of the wedge tip in the 
initial plane. The nominal free stream Mach number at this upstream location 
was set at 2.94, the tunnel total pressure was set at 88.9 psia and the tunnel 
free stream total temperature was 445.0* Rankine. This gave a free stream 
average operating Reynolds number of 1.28E6, based on a refenence length of 
0.667 feet. The tunnel sidewall turbulent boundf'y layer in the initial plane 
was sized to match the measured profiles at the reference survey plane at 
X « 3.5 Inches using a Maise-McOonald initial boundary layer profile 
(ref 11). This gave a boundary layer thickness of 0.55 inches upstream of 
the interaction which was assumed to be uniform across the entire tunnel 
sidewall . 


Boundary Conditions and Turbulence Model 

The PEPSIS computer analysis used for this study has the provision to in- 
dependently prescibe separate boundary condition on each of the bounding sur- 
faces of the computational domain represented in figure 3. In each of the 
four test cases reported in this paper, the McDonald-Camarata turbulence model 
(ref. 11) was applied to all waP surfaces and wall function boundary condi- 
tions were prescribed on the ramp surface downstream of the ramp tip. Upstream 
of the ramp tip, the streamwise component of velocity was held fixed. In the 
region of the ramp tip, the wall function boundary conditions were applied 
using an obstacle logic procedure developed by Buggeln (ref. 6) in order to 
minimize the computational problems associated with the leading edge ramp tip 
singularity. Along the tunnel sidewall, two types of wall boundary condUions 
were studied. The baseline Wall Function case applied wall function boundary 
conditions to each of the bounding surfaces. In addition to this Wall 
Function case, two additional case were studied using no slip boundary condi- 
tions along the tunnel sidewall. The Coarse Mesh case used the identical 
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computational mesh as the Wall Function baseline case where the yplus of the 
first grid point off the tunnel sidewall was 271. This placed the first mesh 
point off tne tunnel sidewall outside the laminar sublayer. A Medium Mesh 
case was also studied In which the mesh points were packed normal to the 
tunnel sidewall such that the first mesh point away from this wall was at a 
yplus of 71, which Is marginally outside the laminar sublayer. To maintain 
the same resolution In the outer region of the flow field, the number of mesh 
points normal to the sidewall was Increased from 40 to 60. Since only half 
the tunnel was considered In this analysis, symmetry boundary conditions were 
applied on the surface passing through the tunnel centerline. To save mesh 
points and consequently to reduce the overall computer run time, symnetry 
boundary conditions were also applied on the surface represented by the tunnel 
celling. This boundary condition was later changed to a wall function 
boundary condition I.e. Tunnel test case, to study the Influence of the 
boundary layer developing along the tunnel celling and to examine the effect 
of the Impinging glancing sidewall shock wave Interaction. To perform this 
calculation, the number of mesh points normal to the ramp surface was In- 
creased from 40 to 80 to resolve the tunnel celling boundary layer. 

The PEPSIS computer design code also has the Input option of Implementing 
either constant enthalpystagnatlon flow conditions or conq)ut1ng the static and 
total temperature field by using the full energy equation. For the series of 
calculations to be presented, constant stagnation enthalpy conditions were 
employed. 
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Con^utatlonal Mesh and Run Time 

A summary of the test cases showing the rr«sh and computing time on both 
the IBM 370/3033 and CRAY I Is presented In table I. The two computing times 
listed under the IBM 370/3033 computer represent two operating modes of the 
NASA Lewis Cohiputlng System. With an Internal Symbol Olctonary (ISO), com- 
putations can be performed In an Interactive mode. This allows the user to 
directly Interact with the computation of the problem while the code Is In 
core. When this option Is removed, the computing time Is reduced by about 40 
percent for the PEPSIS computer code. While the test cases listed In table I 
were not actually run on the CRAY I computer, the times listed are estimated 
from PEPSIS cases that have been run. The CRAY I computing time listed does 
not represent a fully optimized compter code for this machine, and much work 
Is planned In this area. In general, the computing times shown In table I are 
highly competitive In this stage of development of the PEPSIS computer design 
code. 


Calculated Results and Comparison with Experimental Data 

Pitot pressure - The predicted pitot pressure distribution at the 
reference survey plane X - 3.5 Inches Is compared with the measured data for 
the five rake positions In figures 4 and 5. Figure 4 presents a comparison 
between the Wall Function and Coarse Mesh test cases while figure 5 shows the 
effects of mesh resolution by comparing the Coarse and Medium Mesh test 
cases. As can be seen In both figures 4 and 5, a discrepancy between data and 
the analytical prediction can occur because of Insufficent mesh resolution In 
the near wall region when no-sllp conditions are applied. Both the wall 
function results, figure 4, and the medium mesh no-sllp results, figure 5, 
provide very good agreement with the pitot pressure data upstream of the side- 
wall Interaction. While the Initial data applied was sized to match the 
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boundary layer thickness at this reference plane. It becomes quite clear that 
even In a region with no substantial pressure gradient, discrepancies an arise 
with no slip with Insufficient near wall resolution. Figures 5, 6 and 7 
present the same study but In a region where large transverse static pressure 
gradients exist as a result of the Interaction of the shock wave with the 
tunnel sidewall boundary layer. The predictions of pitot pressure when wall 
function boundary conditions were applied at the tunnel sidewall show 
remarkable agreement with data, figure 5. However, even with moderate near 
wall resolution as represented by the Medium Mesh test case, a discrepancy 
occured In the near corner region between the ramp and tunnel sidewall 
surfaces. The differences between the analytical prediction and experimental 
data may be due to any of several possibilities, the most likely being 
Insufficient mesh resolution In the near corner region. This suggests that to 
properly resolve the global flow field In this Interaction, the small length 
scales associated with the corner flows may have to be resolved. It should 
also be pointed out that the near wall and near corner regions are very 
difficult to experimentally survey and every effort should be made to obtain 
accurate data. 

Static pressure - The predicted static pressure distribution as obtained 
by using both wall functions and no-sllp boundary conditions are presented In 
figures 9, 10 and 11 and compared with the measured data at the rake survey 
planes located at X > 5.1 and X > 7.6 Inches. In general, excellent agree- 
ment with the static pressure data was achieved for all the types of applied 
boundary conditions considered In this paper. This suggests that the static 
pressure prediction is not as sensitive to Insufficient near wall resolution 
as the total pressure. The largest discrepancy occured at the rake survey 
plane when the shock wave was well outside the boundary layer and In the free 
stream, YG » 2.75 inches, figure 11. This was caused by the lack of mesh 
resolution in the outer portion of the flow field. 

Figure 12 presents the static pressure signatures calculated on the ramp 
surface and tunnel sidewall. Also shown is a schematic diagram of the com- 
plete tunnel configuration showing the location and orientation of the wedge 
and resulting shock wave. The numerical formation of the shock wave in the 
region of the wedge tip can clearly be observed as well as the shock wave on 
the tunnel sidewall surface. Although not explicitly shown in figure 12, an 
adverse pressure gradient exsists in the near wall region of the wedge surface 
and tunnel sidewall which was measured in (ref. 9). Of special interest also 
is the fact that static pressure gradients are established or the ramp surface 
as a result of the corner flow boundary layer development and that the maximum 
static pressure on the ramp surface occurs just away from the tunnel sidewall 
causing low energy fluid to flow in towards the corner region. 

Yaw angle - The yaw angle is defined as the ratio of the velocity in the 
Y6-d1rection which is parallel to the tunnel sidewall, divided by the velocity 
in the X-direction, see figures 2 and 3. This Is a particularly difficult 
parameter to calculate and measure since it represents the ratio of two veloc- 
ities which approach zero as the wall is approached. Therefore, it would be 
expected to be the most sensitive and the most difficult flow variable to 
match with experimental data. 

Shown in figures 13 and 14 are the effects of applying wall function 
boundary conditions along the tunnel sidewall, figure 13, and the effects of 
Increasing mesh resolution in the near wall region, figure 14. Both these 
comparisons were performed at the X - 5.1 inch survey plane location. In 
general, both the wall function and medium mesh calculational results showed 
excellent agreemtr't with the experimental data, although the wall function 
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results are somewhat lower than the measured results. Excellent agreement was 
also obtained at rake survey planes of X - 7.1 Inches, figure 15 and the 
survey plane X > 7.6 Inches, figure 16, for the Medium mesh calculatlonal 
results. Of Interest Is the fact that the calculated and measured yaw angle 
Increased In the downstream direction, figures 14 to 16, and eventually exceed 
the Inviscid shock wave angle, YG - 1.75 through 3.00 Inches In figure 16. It 
Is also quite apparent that Insufficient resolution of the near wall regl^'n 
when no-si Ip boundary conditions are applied caused discrepancies to occur In 
the calculation of the velocity field of this Interaction. 

Flow Field Visualization 

A detailed picture of the flow field In the tunnel sidewall region is 
shown In figure 17. The computed secondary velocity shown are relative to the 
computational mesh shown In figure 3. In the vicinity of the tunnel sidewall, 
a very strong vortex is established which elongates and increases In strength 
In the downstream direction. This causes very low energy fluid to accumulate 
In the corner region with resulting low wall shear stress in this region. 

This probably accounts for the the very high heating rates measured by Cskam, 
Vas and Bogdonoff (ref. 9). The increasing strength of t.^e sidewall vortex is 
also suggested by the experimental data results since the maximum yaw angle 
increases in the downstream direction, as can be seen in figures 14 to 16. 

The very rapid formation of this sidewall vortex is illustrated in figure 18, 
where the the wedge leading edge region was magnified to Illustrate this 
formation. The overturning that occurs in the sidewall boundary layer results 
from an imposition of the main stream static pressure gradient upon the low 
momentum near wall viscous flow. This overturning results in low energy fluid 
being drawn in towards the tunnel sidewall in the lower ramp region. In the 
outer ramp region, the flow rolls over to form the sidewall vortex. The 
surface velocity vectors, i.e., the velocity vector at the first mesh point 
off the wall is presented in figure 19. As the flow develops in the 
downstream direction, the surface velocity vectors tend to align themselves 
with the shock angle. The shear stress signature that results from this 
secondary flow structure is shown in figure 20. The lowest wall shear stress 
occurs in the near corner region of the wedge surface and tunnel sidewall 
intersection. Also shown in figure 20 is a schematic diagram of the tunnel 
configuration showing the orientation of the wedge and the resulting inviscid 
shock wave. 


Complete Tunnel Calculation 

Of particular importance in the design of two dimensional supersonic 
inlets is the effect of the sidewall vortex system when it impinges on the 
tunnel ceiling and how best to control this interaction. To obtain some 
insights into this benchmark interaction, a fourth or Tunnel test case was run 
where the ceiling boundary layer was also resolved. Figure 21 presents the 
results of this calculation in the form of secondary velocity vectors. At 
station 92, the flow is such that the corner just begins to react to the side- 
wall vortex system. The shock wave in the outer portion of the flow field 
shown has not come into view. By station 98, the original shock wave appears 
in the outer region and the corner flow begins the shock wave reflection 
process. Although the static pressure rise is not sufficient to separate the 
flow on the center portion of the tunnel ceiling, the flow begins to sep- 
arate within the corner region causing a strong shock wave to develope. This 
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shock wave propagates Into the outer region In a diagonal direction as can be 
seen between stations 101 through 116. Thus a hl^ly three-dimensional re- 
flected shock wave develops which separates the flow In the corner region. 
Unpublished experimental data obtained in the 1x1 ft. Supersonic Wind Tunnel 
at Lewis Research Center for similar conditions confirms the existence of this 
separated region. 


CONCLUSIONS 

This paper describes a numerical solution to the glancing sidewall oblique 
shock wave turbulent boundary layer Interaction In three-dimensions by a 
supersonic spatial forward marching computer design code PEPSIS. A test case 
was run to investigate the use of the wall function boundary condition approx- 
imation for a highly complex three-dimensional shock boundary layer Inter- 
action. In general « very good results were obtained with the wall function 
condition boundary approximation. Two additional test cases were run to 
examine the question of near wall resolution when no-sllp boundary conditions 
are used. In the Coarse Mesh test case* the first mesh point off the wall was 
at a yplus of 271.0* while for the Medium Mesh case* the yplus was 71.0. The 
yplus for this latter case was marginally close to the laminar sublayer and 
showed substantially better comparlslon with experimental data* particularly 
the yaw angle data. The differences between the Medium Mesh numerical results 
and experimental data occur In the pitot pressure distribution In the near 
corner region of the flow field. While these differences could be due to any 
of several possibilities* the most likely Is Insufficient mesh resolution in 
the near corner region. This suggests that to properly resolve the global 
flow field* the small lengths scales associated with the flow In the near 
corner region may also have to be resolved. The resulting flow field was such 
that a strong sidewall vortex was established which elongated and Increased In 
strength In the downstream direction. This vortex system caused low energy 
fluid to accumulate In the corner region resulting In very low wall shear 
stress. When this sidewall vortex system Impinged on the tunnel celling sur- 
face, the flow In the corner sepa. ated and formed a highly three-dimensional 
shock wave system which propagate out In a diagonal direction. 

In vgeneral, the PEPSIS computer design code gave excellent results com- 
pared with the sidewall shock wave turbulent boundary layer Interaction data 
of Oskam, Vas and Bogdonoff. Although much work Is planned to Improve the 
fluid flow physic., generalize the geometry and decrease the computing time by 
conversion and optimization on the CRAY I computer* It Is self evident that 
the PEPSIS compter code Is practical for three-dimensional supersonic Inlet 
calculations. 
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TABLE I. - TEST CASES. MESH AND COMPUTING TIME 


CASE 

MESH 

COMPUTING TIME (MIN) 



IBM 370/3033 

CRAY 1# 



ISO 

NO ISO 


WALL FUNCTION 

40x40x120 
(192 000) 

53. 1 

62.0 

7.2 

COARSE MESH 
YPLUS-271 

40x40x100 
(160 000) 

44.0 

26.0 

6.0 

MEDIUM MESH 
YPLUS-71 

40x60x90 
(216 000) 

60.0 

36.0 

9.6 

TUNNEL 

1 

80x40x120 
(384 000) 

104.0 

62.0 

14.4 


^ESTIMATED TIME 
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Figure 4 - Effect of umII functions on the pHot pressure distri 
button in the reference survey plane. X • 3.6 in. 



Figure 5. * Effect of mesh resolution on the pHot pressure distri 
button in the reference survey plane. X - 3.6 in. 
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Figure 6. * Effect of ivall functlonr on the pilot pressure dls- I 

tribution In the survey plane X « 5. 1 in. | 



Fl^re 7. • Effect of mesh resolutlor. on the pitot pres- 
sure distribution in the survey plane X * 5. 1 In. 
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Flgurt 10. * tffmt oT mtsh miolutkin on sUHc prtssurf In 
Wit surety pliw X • M in. 
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Figure W. - EffiKt of mtih itiolution on tht stifle prtssurt dis- 
tribution m thi Sill VYV pliiw X « 4b In, 
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Figure U - Effect of mesh resolution on the yew angle dis- 
trSiution in the survey plane X - 5.1 in. 



Figure 15. - Effect of ,?iesh resolution on the yaw angle dis- 
tribution in the survey plane X • 7. 1 in. 
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Fiqure 16. - tff«ctof mesh resolutkni on the yew angle distrN>u 
tlon in the survey plane X • 7.6 in. 



Figure 16i - Continued 
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Figure 18. - Formjtion of the sidewall vortex. 



Flow 

Figure 19. - Surface veloci^ vectors. 
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Figure JO. - Shear stress signature on the tunnel wall. 
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Figure Jl. - Impingmeni ol the sidewall vorte« on the tur-.el 
celling. 



